Localized matter-waves patterns with attractive interaction in rotating potentials 
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We consider a two-dimensional (2D) model of a rotating attractive Bose-Einstein condensate 
(BEC), trapped in an external potential. First, an harmonic potential with the critical strength is 
considered, which generates quasi-solitons at the lowest Landau level (LLL). We describe a family of 
the LLL quasi-solitons using both numerical method and a variational approximation (VA), which 
are in good agreement with each other. We demonstrate that kicking the LLL mode or applying 
a ramp potential sets it in the Larmor (cyclotron) motion, that can also be accurately modeled 
by the VA. Collisions between two such moving modes may be elastic or inelastic, depending on 
their total norm. If an additional confining potential is applied along with the ramp, it creates a 
stationary edge state. Applying a kick to the edge state in the direction of the ramp gives rise to 
a skipping motion in the perpendicular direction. These regimes may be interpreted as the Hall 
effect for the quasi-solitons. Next, we consider the condensate trapped in an axisymmetric quartic 
potential. Three species of localized states and their stability regions are identified, viz., vortices 
with arbitrary topological charge m, "crescents" (mixed- vorticity states), and strongly localized 
center-of-mass (cm.) states, alias quasi-solitons, shifted off the rotation pivot. These results are 
similar to those reported before for the model with a combined quadratic-quartic trap. Stable pairs 
of cm. states set at diametrically opposite points are found too. We present a VA which provides for 
an accurate description of vortices with all values of m, and of the cm. states. We also demonstrate 
that kicking them in the azimuthal direction sets the quasi-solitons in epitrochoidal motion (which 
is also accurately predicted by the VA), collisions between them being elastic. 

PACS numbers; 03.75.Lm, 05.45.Yv 



I. INTRODUCTION 



Formation of vortices is a well-known manifestation of superfluidity in Bose-Einstein condensates (BECs) - in 
particular, in those which form effectively two-dimensional ("pancake") configurations in appropriately designed 
trapping potentials jlj. If a condensate with repulsive interactions between atoms, which is confined in a nearly 
2D layer by a 2D harmonic (quadratic) potential with trapping frequency w, is set in rotation with frequency 17, 
formation of a multi- vortex lattice is observed in the experiment [2]. The stability of such lattices is limited to < uj, 
as otherwise the centrifugal force empties the region in the center of the trap 3]. Close to the instability threshold, 
formation of a metastable state in the form of a giant vortex, with topological charge m ~ 50, was observed Q. 
On the other hand, it was proposed theoretically Q and implemented in the experiment Q that the instability at 
n > uj can be eliminated if the trapping potential is steeper than harmonic, the simplest possibility being to add a 
quartic term to it [in the critical case of — uj, the linearized version of the respective 2D Gross-Pitaevskii equation 
(GPE) is tantamount to the Schrodinger equation for a charged particle in the uniform magnetic field, which gives 
rise to the Landau levels, see Eq. ^ below]. In many theoretical works, it has been demonstrated that self-repulsive 
condensates can form stable vortices with a multiple topological charge, m > 1, in 2D anharmonic traps 

Dynamics of vortices in BEC with attraction between atoms is different ~ in particular, due to the possibility of 
the collapse in self-attractive media Q . The stability of 2D vortices confined by the harmonic potential was studied 
in detail Q. It was demonstrated that, prior to the onset of the collapse, the vortex with m = 1 is destabilized by 
azimuthal perturbations that split it into mobile localized objects resembling fundamental solitons, if the norm of the 
vortex exceeds a certain critical value. Vortices with m > 2 are completely unstable in the same setting. 

It was also predicted that, under the action of the rotation, BEC with the intrinsic attraction can break the 2D 
axial symmetry by self-trapping into quasi-soliton objects, alias "center-of mass" (cm.) states, characterized by an 
offset of the cm. from the rotation pivot [10.1. I t was concluded that the anharmonicity of the trapping potential is 
necessary for the stability of the cm. states Therefore, the theoretical study of rotating attractive condensates 
trapped in quadratic-quartic radial potentials has drawn attention. Phase diagrams of this model were investigated in 
detail, both in the mean-field approximation [i.e., using the GPE and its linearization for small perturbations, in the 
form of the Bogoliubov - de Gennes equations] |12l ] , and by means of a numerical diagonalization of the many-body 
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bosonic Hamiltonian both approaches producing similar results. Three types of stable localized patterns were 
identified in these studies: vortices with topological charge m — 0,1,2,3,...; crescent-shaped states with a broken 
axial symmetry, that may be realized, for instance, as a superposition of vortices with m — 2,3,4; and the cm. 
states shifted from the rotation pivot. The transition between the crescents and cm. states which feature stronger 
localization is gradual. It is also relevant to mention that crescents built similar to those reported in Rcfs. [ll[ and 
[13] can be made stable in a completely different model, viz., a quasi-linear 2D equation with the harmonic trap whose 
strength is proportional to the total norm of the configuration [14l| (the so-called "accessible-soliton" model, which is 
relevant to nonlinear optics 15]). 

The previous analysis of these states was carried out, in the framework of the GPE, in a numerical form only. One 
of purposes of the present work is to demonstrate that the entire family of vortices, with all values of m, and the 
well-localized cm. states can be predicted, in an accurate form, by a simple variational approximation (VA) (this 
method was first applied to BEG in Ref. [1^]; for a general review of variational techniques for solitons, see Ref. [17|). 
In addition, we find stable pairs of cm. states set at diametrically opposite points. Another objective of this work is 
to consider the motion of cm. states and collisions between them (which turn out to be elastic). The motion is also 
accurately described by the VA. To focus on effects of anharmonic traps, in that part of the work we consider the 
GPE with the quartic radial potential only, which may be implemented in the experiment Q. 

Another setting considered in this work is the above-mentioned critical case of the quadratic axisymmetric trap 
with w = f2. In this case, we demonstrate that quasi-solitons can be found in the lowest Landau level (LLL). Both 
quiescent and moving LLL modes are very accurately described by an appropriately modified VA. The motion, of the 
Larmor (cyclotron) type, is imposed on them by the application of a kick, or by the action of a ID ramp potential. 
In addition, we consider the situation when the ramp acts in a combination with a ID quartic potential. In the latter 
case, we find edge states of the nonlinear LLL modes, and also study their motion induced by the kick, which may be 
interpreted as an effective Hall effect for the quasi-solitons. 

The paper is structured as follows. The underlying two-dimensional GPE and a brief description of the numerical 
method employed to look for stationary solutions (which is based on the integration in imaginary time) are presented 
in Section II. In Section III, we consider the quasi-solitons of the LLL type in the critical model with the quadratic 
axisymmetric confining potential, while the model with the quartic potential is dealt with in Section IV. The paper 
is concluded by Section V. 



II. THE MODEL AND NUMERICAL METHODS 



The mean-field approximation for the 2D condensate trapped in potential U{x,y) and rotating at angular velocity 
n is based on the GPE for the single-atom wave function, t/jix, y, t). In the scaled form, the equation, written in 
rotating coordinates x and y, takes the well-known form Il0l]-|12l]: 

dx"^ 



dt 



V'. 



(1) 



Here, g = —AiTAfas/ {Ma±) > is the effective self-attraction coefficient, with J\f the total number of atoms, 
Us < the scattering length of the attractive interatomic interactions, aj_ the transverse-trapping length, and 
M = J J lijj {x,y) {"^dxdy the norm of the 2D wave function, which is a dynamical invariant of Eq. ([1]). The orbital- 
momentum operator is = i{ydx — xdy). In addition to M, Eq. ([1]) conserves the energy. 



E = 



^]vv]' - f ]^]' + u\^\^ - nrl^ij ] dr, 



(2) 



and the total angular momentum, L — J ip* Lzipdr, if the potential is axisymmetric, U ~ U{r), with 



Steady-state solutions to Eq. ([T|) are looked for in the ordinary form, il){x,y,t) 



<l){x,y), where /i is the 



real chemical potential, while stationary wave function (j){x,y) remains complex. To find configurations realizing a 
minimum of the energy, we used a modification of the known numerical method based on the integration of the GPE 
in imaginary time [To ] (in Ref. [l9| . a similar method was used to generate vortex lattices in a 2D model with the 
repulsive nonlinearity and harmonic confining potential). To this end, we substitute real time t in Eq. ([T]) by — ir, 
introduce additional real variable M(r), and replace Eq. ^ by the following system of the Ginzburg-Landau (GL) 
type: 



dr 
dM 
~d^ 



+g\c^\^ -U{x,y) + nL,+^i [m - M 



72 (Mo - M) . 



(3) 



(4) 
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Here, M is the same 2D norm as defined above, but it is not a dynamical invariant of the GL equations, and Mq 
is the target constant value of the norm in the stationary state to be found, while 71 and 72 are auxiliary positive 
constants. 

Equation @ acts as a negative feedback, which provides for the relaxation of the variable norm, M(r) , to Afo- 
Obviously, stationary states, into which solutions to Eqs. ([3]) and ([4]) relax at t — > 00, also yield stationary solutions 

to GPE ll]), with chemical potential /i = 71 M (r 00) — Mq . 

Coupled equations ^ and ([4]) can be presented in the gradient form, as 



SE 



dr 



dM _ 72 dE 
dr 71 dM ' 



(5) 



where E ^ E - (l/2)7i(Af - Mq)^ + (1/2)71 (A/ - A/)^. This representation demonstrates that Eqs. © and (g]) may 
be regarded as equations generated by the minimization of functional E, under the constraint (Lagrangian condition) 
that the total norm is fixed, M = Mq. Fast convergence of numerical solutions of Eqs. ([3]) and (j4|) was achieved, 
for instance, with the choice of auxiliary parameters 71 = 3 and 72 = 5. On the other hand, if the convergence is 
achieved, the results do not depend on the choice of 71 and 72, which was verified by varying these constant in a 
broad range. The numerical integration was performed by means of the split-step 2D-Fourier method with 256 x 256 
modes. 



III. LOWEST-LANDAU-LEVEL STATES 



A. Quiescent modes 



In the critical case when the confining potential is harmonic, with the respective frequency equal to the rotation 
velocity, U = (1/2) fir^, Eq. ([T]) and respective energy ([2]) reduce to 



1 f.d_ 

2 Vdx 



dy 



(6) 



E 



i|(V-^A)^p-||^|4)dr, 



(7) 



where A — (—Qy,flx). As said above, this system with g — is equivalent to the 2D Schrodinger equation for 
a charged particle in uniform magnetic field O directed perpendicular to plane {x,y). The wave function of the 
corresponding ground state, i.e., the lowest Landau level (LLL), is 



^Jj^A e-'"* exp [- {n/2) {{x - x^f + {y - y^f)] 



(8) 



where xq and j/o determine an arbitrary central position of the particle, and A is an arbitrary amplitude. 

For 5 > 0, a localized (quasi-soliton) solution to nonlinear equation ([6]) may be approximated as an ansatz suggested 
by exact solution ([8]) for the linear equation. 



^ = A e-*^* exp [~a {{x - xof + {y - yof)] , 



(9) 



where A, a and fi are treated as variational parameters. Earlier, an ansatz using a product of the LLL wave function 
and an appropriate function in the perpendicular direction was used in Ref. ^2Q] for the description of 3D vortices in 
rotating confined BECs. 
The norm of ansatz ^ is 

M ^ A^TT/{2a). (10) 
Using this expression, we eliminate A in favor of A/, and then calculate energy ([7]) corresponding to the ansatz: 

gMa 



E = M 



1 Qp- 
a + -il\xl + yl) + — 



2tt 



(11) 
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FIG. 1: (a) A typical shape of the nonhnear localized mode supported by the LLL (lowest Landau level) for M = 5. (b) 
Comparison of the numerically found inverse- width parameter a with the approximation provided by the variational method. 
Other parameters in this figure are g — 1 and fl = 0.5. 



which does not contain /i. To predict the inverse-width parameter a of the locahzed state, we minimize the energy 
with respect to a, by setting dE/da = 0. This yields 

n 1 



2 Vi-gM/(2^)' 

Then, the substitution of this resuh into Eq. (jlOp yields the respective expression for the amplitude. 



(12) 



.9 , , 

A = , 13 

- 9MI(2t:) 

Note that expression ([13]) diverges at gM — 2n, which implies the collapse in the 2D setting due to the self-attraction. 
This collapse threshold is a well-known prediction of the VA 21], which does not depend on the presence of the 
external potential or rotation. 

Figure 1(a) displays a numerically found profile of the central cross section of ip{x,y) for g — 1 and M = 5. The 
dashed curve in the same figure is the Gaussian fitting to the numerical profile, ipatix, y) — Af^t exp [— Q!fit(x — i/2)^] , 
with Afit — 1.62 and agt = 0.824. Note that the central part of the 2D mode is slightly narrower than the Gaussian, 
because M is rather large, inducing self-compression of the mode. The solid curve in Fig. 1(b) displays the prediction 
of the VA for inverse- width parameter a, given by Eq. (|12p . while the chain of circles represent numerical values of 

the same parameter, which were found from an integral expression, M/ 2 J J {{x — xq)'^ + (y — yo)^)|'0(2;, y)\'^dxdy . 

Indeed, if ansatz ([9|) is substituted into this expression, it will yield exactly a. 

Modes with the intrinsic angular momentum, i.e., localized vortices, can also be constructed in this setting. However, 
a detailed analysis of the vortex dynamics is beyond the scope of the present work. 



B. The Larmor motion 



If the localized state placed at (xq, yo) = (0, 0) is kicked by lending it wavenumber ky in the y-direction, ip{x, y) — s- 
ip{x,y) exp^ikyy), the quasi-soliton exhibits rotating Larmor (cyclotron) motion, as shown in Fig. 2(a). If a ramp 
(constant external force) is applied in the x-direction, by adding the extra potential term, C/extra(a;) — —Fx, to 
energy ([7]), the quasi-soliton gets engaged in a drift motion, as shown in Fig. 2(b). These scenarios of the Larmor 
motion of the quasi-soliton are qualitatively the same as exhibited by the LLL localized state in the limit of the linear 
Schrodinger equation {g — 0), as shown in Fig. 2(c). One can check, by means of direct simulations, that a localized 
vortex with vorticity m = \ also exhibits the Larmor motion in the same linear Schrodinger equation. 

The Larmor motion induced by the additional potential can also be explained via the variational method. The 
Lagrangian corresponding to Eq. ([6]), with the addition of the extra potential, is 

L = j (^^(^tV*-V':'/')-^l(V-zA)Vp + ||^|^-f/extra(a:,2;)|^p)dr. (14) 

If ^{x,y,t) is approximated by a generalization of ansatz ([9]), viz., 

i; = Ae-'^^'exp[-a{{x-xo{t))' + {y-yo{t))')] 

X exp[ipx{x - xo{t)) + ipy{y - yo(i))], (15) 
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(a) (b) 




FIG. 2; (a) The Larmor motion of the cm. of the LLL quasi-soliton initiated by the application of kick fey = 2 in the y-direction, 
for g — 1 and M = 5. (b) The drift of the cm. in the ramp potential, t/cxtra = —O.lx, again for g = 1 and M = 5. (c) The 
drift of the cm. under the action of the same tilted potential, but for for g — (in the linear Schrodinger equation). 



the corresponding effective Lagrangian takes the following form, after a straightforward calculation: 



^ = Pxio+Pyi/o- a-^{xl + yl) -^{pI+pI) 



-j-^^ + ^^-^ - ^{PxVa -PyXo) - Ucs {xo,yo) , (16) 
4a ZTT 



where the overdot stands for the time derivative, and 

Uc«{xo,yo,a) = ^ J J Uc^trii{x,y) exp [-a {{x - xo(t)f + {y - yo{t)f)] dxdy. (17) 
The system of the Euler-Lagrange variational equations is then derived from the effective Lagrangian: 



dpx 
dt 


= ^Py 


dU,s 
dxo 


dpy 

~dt 


= -Qp 


X n J 

oyo 


dxo 
~dt 


= Px + 


^yo, 


dyo 
dt 


= Py- 


flxo, 


-1 


gM 




+ t 





(18) 



da 



= 0. 



In the simplest case of the uniform ramp, corresponding to C/cxtra = —Fx, i.e., dUcs/da = 0, the last equation in 
system (|18|) yields the same constant expression for a as given above by Eq. (fT2|) . With constant a (i.e., constant 
width of the LLL mode), the remaining part of system amounts to coupled equations of motion of the second 
order, 

d'^xp _ _ ^^eff 

dt"^ dt dxo 

(19) 

dt'^ dt dya 

With Ucff = 0, this system describes the Larmor (cyclotron) motion, and for UcS = —Fxq, Eqs. (flQl) predict an 
overlap of the Larmor rotation and drift in the y-direction. The drift motion, in the direction perpendicular to the 
ramp, is an analog of the ordinary Hall effect in solid-state physics (similarities between the soliton dynamics in BEC 
and the quantum Hall effect were discussed in various contexts, see, e.g., Ref. [l^ and references therein). Additional 
analogies to the Hall effect are considered below. 

If two LLL quasi-solitons are kicked in opposite directions, by imparting wavenumbers ±fc,y to them, the Larmor 
motion of the solitons eventually results in a collision. If the norm of each quasi-soliton is not too large, the collision 
seems completely elastic, and the two objects keep moving in closed trajectories which together form a "figure of 




FIG. 3: (a) Trajectories of points of maxima of |?/>(a::,y)| in half-planes x > L/2 and x < L/2 for a pair of LLL quasi-solitons 
with norm M = 1, kicked by ky = ±2. The solitons survive multiple elastic solutions. The arrows shows the initial positions 
and velocities of the two solitons. (b) The same for M = 5. In this case, the collision results in the fusion and collapse of the 
LLL quasi-solitons. 



eight", surviving multiple collisions, as shown in Fig. 3(a) for M = 1, g = 1 and ky = 2. On the other hand, if 
the norm of each LLL mode is larger, the collision is inelastic, reducing the wavenumbers of the quasi-solitons, and 
thus making the radius of the Larmor motion smaller. This case is illustrated by Fig. 3(b) for M = 5. Because the 
respective total norm, Mtot — 10, exceeds the collapse threshold, the two quasi-solitons merge into a single collapsing 
object, after several consecutive collisions. Inelastic collisions between 2D solitons which are too "heavy" and can 
also suffer the merger and collapse were reported in different settings, such as quasi-lD guiding channels [23|. 



C. The Hall effect and edge states for the quasi-solitons 

As mentioned above, the analog of the Hall effect for matter- wave solitons is a subject of considerable interest [2^ . 
To study it in the present context, i.e., as a matter of fact, to consider the corresponding edge states of the LLL quasi- 
solitons, we combine the ramp with a weak holding quartic potential in the a;-direction, by taking f/oxtra = —Fx + bx'^, 
with small b > 0. 

The VA can be used in this setting too. To this end, we approximate ip by an anisotropic Gaussian ansatz, 

^ = A e-'^'* exp[-a{x - xof - (3{y ~ yof] 

exp[ipx{x - XQ{t)) + ipy{y - yoit))], (20) 

cf. expression ([T5|) . The substitution of ansatz ^U\j in Lagrangian ([H]) yields 

^ = PxXo + Pyi/o - ^{a + (3) - ^{xl + yl) ~ ^{pI + pI) 

(a ^ ^) ^ ^(P^yo -Pya;o) - Ucs {xo,yo) , (21) 

where, this time, the norm is M — A^tt/ (2^/a/3) , and 

C/off {xo,yo) = ~Fxo + b{xl + yl f + ^{3x1 + y^) + ^{xl + iyl). (22) 
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FIG. 4; (a) The contour plot of Re{i>{x,y)} of the edge state, for 5 = 1, M = 5, F = 0.5 and b = 0.0005. (b) The profile of 
\tp{x,y)\ at y — L/2. (c) Ke{tp{x,y)} at a; = 21.2. The two latter panels include comparison with profiles predicted by the 
variational approximation. 



cf. Eqs. (|16p and p?)) . The respective equations of motion are written as 



d^yo 
dt^ 



2n 



dyo dUcs 



gM 
in 
gM_ 
47r 



dt 

dxo 



dxo 



dt 

a 
fa 



dyo ' 






at/off 




da 




dU,s 


8/32 





(23) 



0, 
0. 



Stationary solutions to Eqs. ((23| are determined by algebraic relations, 

yo = 0, Abxl + 5(3/a + l/^)xo = F 



(24) 



- (1/2) + 5Af/(47r) V/V" + f^V(8a^) + 36x^/(2^^) = 0, (25) 

- (1/2) + gM/(4^) f^V(8/3') + foa;2/(2/32) = 0. (26) 

In the linear limit (3 = 0), these relations yield a = ^JVfi /A + 36xq and /? = ■\/il2/4 + 5xq, with wavenumber in the 
y-direction being py ~ Q,xq > 0. 

As an example, we refer to the LLL mode in the edge state, that was found in the numerical form for M = 5, 
F = 0.5, b = 0.0005, (? = 1, and fl = 0.5. Fitting this mode to Gaussian ansatz (PO)) . the corresponding parameters 
were numerically evaluated as xq = 6.225, a = 0.724, /3 = 0.679, A = {Aa(3 /ttY/'^Vm = 1.99, and Py = Q.xq. 
Figures 4(a) and 4(b) display, respectively, the contour plot of Re{-0(a;,y)} and cross-section profile |V'(a;,2;)| at 
y — which is compared to its counterpart predicted by the VA through Eqs. ([20]) and ((25)) . In addition. Fig. 
4(c) displays Re{^(x,?/)} at a; = 21.2, which is compared to the respective approximation provided by ansatz (|20|). 
Aexp{-/3(y - L/2)2} cos{pyy). 

If the edge state is kicked with wavenumber ky, the quasi-soliton exhibits drift motion. For parameters identical 
to those in Fig. 4, in Fig. 5 we display the trajectory of its cm., initiated by the kick with ky — 1, and the 
same trajectory as predicted by variational equations (|23p . These results demonstrate that the VA provides quite a 
reasonable accuracy for the description of dynamical states, as well as for static ones. 



IV. THE AXISYMMETRIC QUARTIC POTENTIAL 



A. Quiescent states: vortices, quasi-solitons, and crescents 



In the rest of the paper, we focus on the situation opposite to that considered in the previous section, namely, the 
setting with the quartic axisymmetric potential, U{x,y) = hr^ , where 6 is a small positive constant. The purpose 
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FIG. 5: The trajectory of the drift motion of an edge LLL mode, kicked with ky = 1, for g = 1, F = 0.5, b — 0.0005, and 
M = 5, as obtained from direct numerical simulations of the GPE (the solid curve), and the counterpart of the same trajectory, 
predicted by the variational approximation (the dashed curve). 



of the consideration of this model is to highhght the dynamics on vortices and quasi-solitons under the anharmonic 
confinement. The respective variant of the GPE is [cf. Eq. ©] 



1 



(27) 



where is the 2D Laplacian, and the angular-momentum operator is — i {xdy — ydx)- 

Similar to Eqs. ([3]) and (|4]), steady-state solutions to Eq. ([T]) are looked for through the simulations of relaxation 
in the following modified system of the GL type: 



dr 



dr 



(28) 
(29) 



Here, M is the same 2D norm as defined above, and Mq is the target constant value of the norm in the stationary 
state to be found. 

Using the remaining scaling invariance of Eq. ([T]), we fix normalizations by choosing h — 0.002 (a small value of 
trapping coefficient h is necessary to allow the condensate enough room to evolve), and Afg — 4. Axisymmetric vortex 
solutions to Eqs. (p8| and (|29p with topological charge m are sought for as </> = r^e^™^ Rm{r,T), where the reduced 
amplitude function, i?^, satisfies the following equations: 



dRm 



2m+l d 
r dr 



br^ Rm + .9r2"i?f„ + mURra + a{M - M)R,n, 



(30) 



— = i3(Mo - M), Af(T) = 27r / r^"'+'Rldr. (31) 
ar Jo 

Once stationary solutions were found, their stability was examined by means of computation of the respective eigen- 
values, using the Bogoliubov - dc Gennes equations, i.e., the linearization of Eq. ([T]), for small perturbations built 
(for given m) as superpositions of components with vorticities mil (such perturbation modes turn out to be most 
dangerous in terms of the instability). 

The results, in the form of a stability diagram in the plane of (Q, g) for m = 0, 1, 2, 3, 4, 5, are presented in Fig. 6(a). 
Each stability domain for m > 1 is bounded by two curves, which are generated by critical perturbation eigenmodes 
that are found to be, respectively, real and imaginary, with respect to unperturbed amplitude functions Rmir). This 
stability diagram is qualitatively similar to those reported in the model with a mixed quadratic-quartic radial trap 
[l2,[i3l- For given > 0, vortices with m < can be found too, but they all are unstable. This instability can be 
readily explained by the fact that the Coriolis term in the vortex' energy, which is proportional to — mf2, is positive 
for m < 0, see Eq. ([55]) below. 

While the increase of fl at fixed self-attraction coefficient g leads to the transition to vortices with larger m, the 
vortex states become unstable with the increase of g, being replaced by crescent-shaped ones, as shown in Fig. [71 In 
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FIG. 6; (a) The numerically generated stability diagram for vortices with different values of topological charge m, in the plane 
of the rotation frequency and strength of the attractive interaction, Q and g, in the model with the quartic confining potential, 
(b) Regions where the variational approximation predicts that vortices with charge m or cm. states, alias quasi-solitons (the 
black area), provide for a minimum of the energy. 
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FIG. 7: Generic examples of stable matter-wave patterns in the rotating condensate confined by the quartic potential are 
shown by means of contour plots of \Kc{(f>{x,y)}\: (a) a vortex with m = 4, for g = 0.01, — 0.35; (b) a crescent, for 
g — 0.001, Q = 0.5; (c) a strongly compressed cm. state (quasi-soliton) , for g — 0.5, Q = 0.5. 



Refs. [T2| and similar patterns were interpreted (in models with the local and nonlocal nonlinearity, respectively) 
as superpositions of vortices with three different values of the topological charge, viz., given m and m ± 1. This fact 
agrees with the above-mentioned finding that perturbation modes which can destabilize a given vortex carry vorticities 
mil. Further increase of g leads to a reduction of the crescent's length and its gradual compression into a strongly 
localized cm. state, alias quasi-soliton, which carries an intrinsic phase gradient along the azimuthal direction, see 
Fig. Ulc); note that the value of g corresponding to the quasi-soliton falls below the collapse threshold. The shift of 
the cm. state from the rotation pivot increases with fl. This quasi-soliton is similar to the edge state of the LLL 
type shown above in Fig. 4. 

In the limit of large m, the stationary version of Eqs. (I30p and pip give rise to a simple asymptotically exact 
solution, which does not depend on g and b: 



Rmir) = ^/Mon"'+'^ I (7rm!) exp {~^r^ j^) 



(32) 



(recall we here fix the norm as Mq = 4). This solution, and numerical results obtained for finite to, suggest to approxi- 
mate the stationary solution for vortices by ansatz = ^gime^m g-j^-p (^_q,j,2'^ ^ whose norm is Mq — nml (2a)~^™^^-' A^. 
The substitution of the ansatz in Eq. ([2]) yields the corresponding expression for the energy, 



E 



(TO + 2)(m + l)6 
4q!2 



+ a{m + 1) - {2my.Moga/ 



■)2jn+l 



(to!) TT 



(33) 



Then, width parameter a for the solution sought for is determined by the minimization of the energy, dE/da = 0, 
which yields 



b{m + 2)(m + 1) 



m + 1 



{2m)lMog 

22"+l(TO!)27r 



(34) 



On the other hand, a quasi-soliton with the cm. located at distance xq from the pivot, see Fig. [TJ^c), may be 
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FIG. 8: Parameters of stable quasi-solitons trapped in the quartic axisymmetric potential, as found from the numerical 
computations (circles) and predicted by the variational approximation (curves) for g = 1: (a) the shift from the rotation 
pivot; (b) amplitude; (c) intrinsic wavenumber. 



approximated by the anisotropic ansatz, which is similar to the one used above in Eq. (|20p . 

= Aexp [iky — (^a{x — xq)'^ + , (35) 

with norm AIo = ttA^/ (2^/a/3^ . If this cm. state was generated by an mistable vortex with charge m, comparison of 
the azimuthal phase gradients suggests that k = m/xo. When substituted in Eq. fj]), ansatz ([35]) yields 

E a + l3 + k^ ^, gV^Mo b \3/3^ + 2a(3{l + 12(3x1) + + 8f3xl + 16(3^x1,)] , ^ 

Mo^ 2 l^fcxo--^^ + . (36) 

Values of the variational parameters are predicted by equations dE/dx^ = dE/dk ~ dE/da = dE/dj3 = 0, which 
yield, in particular, 

k = nxa, xl = [n'^af3 - 5(3/3 + a)] / {Abaf3) , (37) 
if xq 7^ 0. Another solution, with xq — (an isotropic soliton sitting at the center, which is stable for small Q) has 

1 /3 

k = and a = /3 = [{2nb)/{2n — Afo)] . Further analysis of the VA solutions (without fixing Mq = 4) demonstrates 
that they predict the collapse (nonexistence of solutions) at Afo > 27r, in accordance with the known variational result 

Using the VA solutions and expressions ((33)) and ((36|) . we have identified, as shown in Fig. [D^b), regions in parameter 

plane (^^,5) where vortices with particular integer values of m, or the cm. state provide for the minimum of the 

energy, i.e., determine the ground state. Comparison with the numerically identified stability regions for vortices in 

Fig. El^a) demonstrates that the VA predicts transitions between different values of m with the increase of quite 

accurately. The discrepancy in Fig. [6] between the numerical and variational plots in the direction of g has an obvious 

reason: the VA does not take into regard the other species of the localized states, viz., crescents, alias mixed- vorticity 

states [see Fig. El^b)]. In fact, crescents have their own domain of the energy dominance, between those of the vortices 

and cm. states. Note that solution ([32]) . which is asymptotically exact for m ^ 00, and expression ([33)l predict the 

equality between energies of vortices with m and m + 1, i.e., borders between the energy-dominance areas of these 

1/3 

states (for g ^ 0), at ~ (46m) . To directly test the accuracy of the VA, in Fig. 8 we display a comparison 
between characteristics of the quasi-solitons, viz., cm. offset xq, amplitude and wavenumber fc, as found from 
numerical results and predicted by the VA. 

The variational results presented above suggest that the model may also support a pair of cm. states placed at 
diametrically opposite points. Indeed, Eq. ([57)1 gives rise to two roots, = ±^\p?-af3 — 5(3/3 -f a)\ / (46a/3), which 
correspond to opposite values of y-wavenumber k. Such stable quasi-soliton pairs can be readily found from the 
numerical solution, see a typical example in Fig. [51 In the experiment, the pair can be created, for instance, by 
originally adding a strong blue-detuned (repulsive) light sheet which cuts the circular trap into semi-circles. After 
two cm. states have self-trapped, the sheet may be turned off, to restore the axial symmetry of the trap. 



V. THE MOTION AND COLLISIONS OF QUASI-SOLITONS 



Although they were obtained above as quiescent solutions, the localized cm. states can be readily set in motion by 
the application of tangential kick exp(iqy) with wavenumber q, similar to how this was done above for the LLL localized 
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FIG. 9: A generic example of a stable pair of two strongly compressed cm. states found for g — 0.5, Q = 0.5. 

(a) (b) 




FIG. 10: (a) The trajectory of the cm. motion of a quasi-soliton kicked with wavenumber g = 1, for g — 1, fl — 0.5, in the 
model with the quartic axisymmetric trap, (b) Its counterpart predicted by the variational approximation. 

modes. As a result, the quasi-soliton exhibits rotary motion, following a trajectory in the form of an epitrochoid, see 
a typical example in Fig. [TOl 

Following the approach elaborated above for the moving LLL modes, we present a similar analytical description of 
the kicked quasi-soliton in the model with the quartic confinement. To this end, we adopt the same ansatz (jisp as 
used above, and use the Lagrangian corresponding to Eq. ([27|) . 

L = j (^-{M* - rti') - ^IV^P + nri.^ + fl^l" - \hAM''^ dr, (38) 

cf. Eq. p^ . The substitution of the ansatz in this Lagrangian and straightforward calculations lead to the following 
equations of motion for the cm. mode: 

dt^ dt dxo ' 

(39) 

d^yo _ _2n— - ^^"'^ 

dt'^ dt dyo ' 

with Ucff{r) = — (l/2)f7^r^ + br"^ + (26/a)r^. Figure [TUf b) displays a counterpart of the numerically found trajectory 
from panel [TUTa) . as produced by Eqs. (|39p for g ~ 1,6 = 0.002, and a — 0.532. The initial velocity is dx^/dt = 
0, dyo/dt = 1. It is seen that the VA matches the numerical findings very well in this case too. 

The possibility of the motion of the cm. states suggests to consider collisions between them, also in analogy 
with what was done above for the LLL modes. In Fig. [TI] we display a typical example of the collision, which is 
generated by applying kicks with q = ±0.5 to identical quasi-solitons with their cm. placed, originally, at diametrically 
opposite points (in the experiment, such an initial configuration can be created as outlined at the end of the previous 
subsection). The kicked solitons are no longer identical because term fiL^ in Eq. ((27| breaks the symmetry between 
the clockwise and counter-clockwise directions of the rotation. Nevertheless, the collisions are elastic, with the quasi- 
solitons recovering their shapes after the collision. 
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FIG. 11: Snapshots of |Re {^'(a;, j/, t} |, at t = 10 (a), t = 50 (b), and t — 90 (c) illustrating an elastic collision between 
counter-rotating cm. states kicked by g = ±0.5, at <? = 0.5, f2 = 0.5. 



VI. CONCLUSION 



We have revisited the 2D model of rotating BEC with attraction between atoms. Two different situations of special 
physical interest were considered: the one with the critical strength of the quadratic confining potential, and the purely 
quartic axisymmetric trap. In the former case, the linear limit of the GPE (Gross-Pitaevskii equation) is tantamount 
to the Schrodinger equation for a charged particle moving in the uniform magnetic field. We have demonstrated that 
the action of the self-focusing nonlinearity on the localized state corresponding to the wave function at the LLL (lowest 
Landau level) gives rise to stable quasi-solitons. These states, both quiescent ones and those set in motion by the 
kick, or under the action of the ramp potential, are very accurately described by the VA (variational approximation). 
We have also considered the situation when an external weak ID quartic potential acts in combination with the ramp, 
which gives rise to edge states emulating the Hall effect in terms of the matter-wave quasi-solitons. 

In the case when the axisymmetric trap is represented by the quartic potential, we have developed the VA which 
provides for an accurate description of two species of stable localized states in the model, namely, vortices with an 
arbitrary value of the topological charge, and cm. modes shifted off the rotation pivot, alias quasi-solitons. Stable 
states in the form of two cm. modes placed at diametrically opposite sites were found too. The other species, 
crescents, was obtained in the numerical form. It was also demonstrated that kicking the cm. state in the tangential 
direction sets it in motion along an epitrochoidal trajectory, and collisions between such solitons are elastic. The 
motion of the kicked quasi-soliton in the latter situation is also accurately predicted by the VA. 
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